Gravitational force in an infinite one-dimensional Poisson distribution 
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We consider the statistical properties of the gravitational field F in an infinite one-dimensional 
homogeneous Poisson distribution of particles, using an exponential cut-off of the pair interaction 
to control and study the divergences which arise. Deriving an exact analytic expression for the 
probability density function (PDF) P{F), we show that it is badly defined in the limit in which 
the well known Holtzmark distribution is obtained in the analogous three-dimensional case. A well 
defined P{F) may, however, be obtained in the infinite range limit by an appropriate renormalization 
of the coupling strength, giving a Gaussian form. Calculating the spatial correlation properties we 
show that this latter procedure has a trivial physical meaning. Finally we calculate the PDF 
and correlation properties of differences of forces (at separate spatial points), which are well defined 
without any renormalization. We explain that the convergence of these quantities is in fact sufficient 
to allow a physically meaningful infinite system limit to be defined for the clustering dynamics from 
Poissonian initial conditions. 
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I. INTRODUCTION 



In attempting to understand better the complex dy- 
namics of classical non-relativistic self-gravitating sys- 
tems relevant in astrophysics and cosmology, it is of in- 
terest to study simplified toy models in one dimension. 
For systems of a finite number of particles, the evident 
such model consists of particles on a line interacting by 
forces independent of their separation (or, equivalently, 
infinite parallel sheets embedded in three dimensions). 
This so-called "sheet model" has been quite extensively 
studied in the literature (see, e.g., [l-Q and references 
therein) for its interest both as a toy model for gravity in 
three dimensions and, more generally, as a toy model for 
systems with long-range interactions. In the study of the 
formation of structure in the universe in cosmology, the 
problem may be well approximated, over a large range of 
length and time scales, by a variant of the simple three- 
dimensional (3D) Newtonian problem: particles belong- 
ing to an infinite distribution of particles evolve under 
equations of motion which are formally identical to those 
for a finite system, up to simple modifications which take 
into account the expansion of the universe (see e.g. [1,Q). 
A few groups of authors EMI have proposed different 
variants of the simple one-dimensional (ID) sheet model 
to mimic these equations. We have underlined in a re- 
cent paper [2^ that a fundamental question about any 
such model, just as in three dimensions, is whether the 
gravitational force term — which is simply the infinite 
sum representing the force exerted by all other particles 
on the given particle — is well defined in the class of 
infinite distributions one wishes to study (which will rep- 
resent the initial conditions for structure formation in 



cosmology). In [20| we have rigorously shown that, with 
an appropriate prescription, the infinite sum in one di- 
mension results in a finite and simple expression for the 
gravitational force acting on particles, in a specific class 
of infinite point distributions: infinite lattices subjected 
to a class of stochastic perturbations. In this paper we 
consider this issue of the definedness of the ID gravita- 
tional force in an infinite and homogeneous Poisson parti- 
cle distribution, and then the related (but not, as we will 
discuss, identical) question of whether an infinite system 
limit for the dynamics of gravitational clustering can be 
usefully defined in this case. 

The approach we use to the question of the definedness 
of the force follows the seminal work of Chandrasekhar 
dating back to 1943 |^]: the infinite point distribution is 
described by a stochastic point process in infinite space, 
and one considers the statistical properties of the grav- 
itational field at an arbitrary point which is itself then 
a stochastic quantity. This approach has been adopted 
and generalized in various cases by ourselves and our col- 
laborators as well as other authors, in the treatment of 
both gravitational and other forces [l^-fl^. In his origi- 
nal work Chandrasekhar derived the result for the PDF 
of the total gravitational field generated by a 3D infi- 
nite homogeneous Poisson distribution, which is given 
by the so-called Holtzmark distribution. It is the general- 
ization of this case to one dimension which is the subject 
of the present paper. We note that this question has 
been treated also in a recent paper [2^. The results we 
present here extend considerably, and allow us to explain 
the physical meaning of, those presented (for one dimen- 
sion) in this latter paper. 

We recall at the outset that to define the gravitational 
force in any dimension in any spatially homogeneous in- 
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finite mass distribution with non-zero mean density it is 
necessary (but not necessarily sufficient) to give a pre- 
scription for its calculation. This is the case because the 
gravitational pair force in d dimensions, which decays as 
r^-d- with the separation r, gives a force on any parti- 
cle due to the non-zero mean density which is badly de- 
fined. In three dimensions the appropriate prescription 
(for the cosmological problem) consists in the so-called 
"Jeans' swindle" (see, e.g.,[33, HJ): the contribution of 
the mean density is removed and only density fluctua- 
tions are taken as the source for the gravitational field. 
As discussed in [s^ this "swindle" is a (well-defined) 
mathematical regularization of the gravitational force 
acting on particles in the infinite system limit (at con- 
stant non-zero mass density), which can be most simply 
re-stated as the prescription that the force on each parti- 
cle be calculated by summing symmetrically about it (e.g. 
in spheres centred on the point). In his calculation of the 
PDF in three dimensions Chandrasekhar adopted, albeit 
without explicit discussion, the prescription of summa- 
tion in such spheres, and other calculations in three di- 
mensions (for other distributions) have done the same. In 
24| one of has shown that the adaption of this prescrip- 
tion to one dimension, i.e. summation in a symmetric in- 
terval about each point, gives a well defined force PDF, in 
a class of infinite perturbed lattices, only for pair forces 
which decrease with separation, but not for the (sepa- 
ration independent) ID gravitational force. In 20| we 
have shown, however, that by adapting a smooth version 
of this prescription — calculating the total gravitational 
field as the limit of the sum for an exponentially screened 
gravitational pair interaction when the screening vanishes 
continuously — we obtain a well defined PDF also for the 
gravitational interaction. Thus while in three dimensions 
the Jeans' swindle can be implemented in various differ- 
ent ways with the same finite results in a wide class of 
stochastic point-mass distributions, the use of a smooth 
version of it (rather than a sharp "top-hat" version) has 
to be preferred in one dimension. This is due to the fact 
that 20] even in a wide class of very uniform stochastic 
point process, stochastically perturbed lattices, a top-hat 
prescription leaves fluctuating and non-converging (i.e. 
undetermined) boundary contributions which are erased 
by a smooth regularization. It is this formulation of the 
Jeans' swindle which we will employ here, although one 
can, as we will discuss, in fact equally use the "top-hat" 
prescription and obtain equivalent results in the specific 
case of a Poisson distribution. 

The first result of this paper is an exact expression 
for the PDF of the screened ID gravitational force in an 
infinite Poisson distribution, derived in Sec. IIIIl More 
specifically we give the exact expression for the cumu- 
lants of the PDF. Using this expression we show that the 
PDF is, as can be anticipated from a simple qualitative 
analysis we give in SecHH badly defined in the limit that 
the inverse screening length /i goes to zero. This in three 
dimensions corresponds to the limit in which instead the 
well-defined Holtzmark distribution has been obtained by 



Chandrasekhar . 

In Sec. IIVI we discuss two other ways in which the in- 
finite range limit may be taken, both involving a renor- 
malization of the other system parameters (mean density 
and coupling strength), and give the (different) asymp- 
totic form of the PDF obtained in each case. The mean- 
ing of these renormalisation procedures is explained in 
Sec. |V] by means of an exact calculation, again for the 
screened interaction, of the force-force correlation func- 
tion. This shows that the non-trivial renormalisation 
procedures actually give a spatially uniform force (ac- 
celeration) field, as in the corresponding limit only the 
(initially) divergent contribution to the force from par- 
ticles infinitely far away survives. This means that the 
dynamics under the renormalised forces corresponds to 
a trivial (albeit stochastic) rigid translation of all parti- 
cles, while the spatially varying component of the forces, 
which would lead to non-trivial relative motions (i.e. in 
this case, clustering) vanishes. 

This discussion leads us naturally to focus on the fact 
that the spatially non-trivial part of the correlation func- 
tion is in fact well defined in the simple (unrenornialized) 
infinite range limit. We show in Sec. IVII that this is re- 
fiected more generally in the fact that the differences of 
forces between points at some fixed distance is manifestly 
well defined and independent of the distribution outside 
the interval they enclose. We derive an exact expres- 
sion for the PDF (specifically, again, for its cumulants) 
of the difference in the (unscreened) gravitational force 
between two points in space, both without and with the 
Jeans' swindle (implemented either with a top-hat or the 
screening prescription). As we are interested in the clus- 
tering dynamics manifested in the finite system, i.e. rel- 
ative motion of particle initially contained in any finite 
region, we discuss whether the definition of the relative 
force PDF is sufficient to make the infinite volume limit 
meaningful. We argue that, if the Jeans' swindle is em- 
ployed, the limit can indeed be defined. 



II. ID GRAVITATIONAL FIELD IN A POISSON 
DISTRIBUTION 

Let us consider a homogeneous Poisson particle distri- 
bution (i.e. a random particle distribution) [26l . Issj on 
the interval [—L,L'] with average density no. In other 
words it is characterized by a microscopic density 



N 

n{x) = 5{x — Xi) 

i=l 



(1) 



where Xi is the position of the i particle, and such that 



{n{x)) = lim 



N 



- no > 



L.L'^oo L + L' 

{n{x)n{x')) = Uq + nQS{x ~ x') , (2) 

where (•) means the usual ensemble average. The sec- 
ond relation says that there is no correlation between 
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the positions of different particles. More generally the 
joint probability density function (PDF) of the positions 
of all the particles is simply 

Pn{xi,...,xn) = (^j-^j-j^ (3) 

The ID version of the gravitational field generated at y 
by a particle at x may be written 

X — y 

fix -y)^ 9] 4 g ■ sgn{x - y) . (4) 

\x - y\ 

The coupling coefficient g gives the intensity of the inter- 
action^. The total field at the point y is therefore given 

by 

F{y) = ■^(^^ -y)= d.^ n{x)f{x - y) . (5) 
i=i •'-^ 

We want to study the statistical properties of this quan- 
tity in the thermodynamic limit N, L, L' — >■ oo with the 
sole constraint N/{L + L') uq. It is clear that in this 
limit the integral Eq. ([S]) is, just as the analogous one in 
3D gravity, ill defined as n{x)f{x — y) is not integrable at 
large |x| in almost any realization of the Poisson particle 
system. 

It is simple to show by the following simple handwav- 
ing argument the result we will find rigorously in the next 
section: while in 3D the usual "Jeans' swindle" suffices 
to make the typical force in a Poisson distribution well 
defined in the infinite system limit, this is not the case in 
one dimension. We recall that Jeans' swindle consists in 
considering as the source of the field not the density field 
n{x), but the density fluctuations 5n(x) = [n{x) — uq] 
from the mean. For any central force, this is equivalent 
to considering the whole n{x) as the source, but summing 
symmetrically with respect to the point y where the field 
is evaluated, e.g. in spheres centered on this point. In 
the introduction of Chandrasekhar's celebrated deriva- 
tion ^21] of the PDF of the total gravitational field in a 
3D homogeneous Poisson particle system, it is this latter 
procedure which is implicitly followed. In one dimension 
this corresponds to taking^ L' = L -I- 2y in Eq. ([5]). To 
see the diff'erence between the two cases (in one and three 
dimensions), it suffices to estimate, once the Jeans' swin- 
dle is adopted, the order of contribution AF^n to Eq. ([5]) 
coming from sources at a distance greater than R from 
the point y. In doing this we approximate the integral 
in Eq. ([S]) with a sum over shells defined by the the se- 
quence of radii i?„ = 2"i? with integer n = 0, 1, 2, i.e., 
radii equally spaced on a logarithmic scale. It is known 



^ It is equal to 27rGS when we derive the model from infinite par- 
allel sheets in three dimensions, where G is Newton's constant 
and E is the surface mass density of the sheets. 

^ This prescription is adopted also in [29II . 



that the typical density fluctuation Sny in a given vol- 
ume V is, for a homogeneous Poisson point process in 
any spatial dimension, of order V^^^'^. In particular the 
volume of the (n + if" shell is K = A(2"i?)'^ where A 
is a geometrical prcfactor depending on d but not on n. 
Therefore we can say that the typical density fluctuation 
in the {n + if" shell is 6n{Rn) ~ [^(2"i?)'^]-i/2. Given 
that the pair force between two particles at distance r is 
of order r~'^^^ in d dimensions, it follows that at suffi- 
ciently large R we can approximate Eq. ([5]) by 

AF>« ^ £ A(2"i?)'^ [^f "^fJ_;'^' ~ C{d)R"^/'+' , 

(6) 

with C a positive constant depending on d. Let us start 
by analyzing Eq. ([6|) for d — Z. We can simply ver- 
ify that C{d) is finite and therefore such a contribution 
to the total force is finite for any R and vanishes for 
i? — > 00. Thus, when one sums in spheres about a given 
point, the typical force on a particle converges. This is 
the fundamental reason why, once the Jeans' swindle in 
this form is adapted, the PDF of the gravitational force 
is well defined (and given, as derived in [2l|, by the so- 
called Holtzmark distribution). On the other hand, in 
d = \ one can see that C{d) — -fcxo, i.e., at any finite R 
the quantity AF^ji is divergent. Moreover the R depen- 
dence is pathologically increasing with R. This means 
that, even if one sums symmetrically, the contribution 
from fluctuations around the mean density due to far 
away regions always give the dominant diverging contri- 
bution to Eq. ([5]). In the following we will demonstrate 
this result more formally, showing that the PDF of the 
total field F calculated in a symmetric window of size L 
becomes broader and broader as L increases, vanishing 
for any finite value of F in the limit L — > (X), i.e., the 
force F is an ill defined and completely undetermined 
stochastic quantity even in this symmetric limit. 

III. PDF OF THE REGULARIZED FORCE 

In order to study in a controlled manner the statis- 
tical properties of the ID gravitational force which, as 
anticipated, may be badly defined, we follow a procedure 
like that adopted often in the context, notably, of quan- 
tum field theory: we introduce a regularization and then 
study the behaviour of relevant physical quantities in 
the limit that this regularization is removed. The Jeans' 
swindle itself can, as discussed in [s^], be considered to 
be such a scheme: in the usual "top-hat" implementa- 
tion the regularization parameter would then be the size 
of the symmetric interval (or radius of the sphere in three 
dimensions) in which one sums. One can equally consider 
a smooth version in which the symmetric sum is imple- 
mented by screening symmetrically the interaction. As 
discussed in the introduction, it is the latter procedure 
we adopt here, as we have shown in [20] that this form is, 
in one dimension, actually preferable to the top-hat form 
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because it gives, in a class of more uniform distributions 
than those considered here, a well defined (and physically 
meaningful) finite force where the top-hat regularization 
does not. We thus consider the pair interaction exerted 
by a particle at x on another at y 



^(x-y)=5[sgn(x-y)]e-'^l^-^l. 



(7) 



introducing a cut-off length characterising an expo- 
nential screening of the "bare" gravitational interaction. 
We will take the limit /j, — 0+ at the end of our calcula- 
tion of physically relevant quantities, and specifically, in 
the next section, to find the existence conditions for the 
PDF of the total force in this limit. 

Given the distribution of particles in the interval 
[— L,i'] defined by Eqs. ([T]), the total field at y is 



N 



L' 



dxn{x)ff,{x - y) . (8) 



Without loss of generality let us now fix y — 0. Using 
Eq. ©, the PDF of = F^{0) given N,L and L' can 
be written as^ 



L N 

.5 



dxj 



U 



N 



F,-Y.U{x,) 



By using the identity 

6{z) 

we can write 



+ 00 



dq 
2^ 



— " X 
2tt 



(9) 
(10) 

(11) 




7 exp 



{-igq[sgn{x)]e ^1^1 } 



N 



In order to take the limit {N, L, L') — >■ oo with N/{L 
L') = no > 0, we start by writing 



dx 



-L'L + L 

j-L' L + L 



- exp|-i5rg[sgn(a;)]e 

'^^ 7 (l -cxp|-ig(7[sgn(a;)]e"^l^l|) . 



Since for /i > we have 

f + 00 



dx {l - exp ^-igq[sgiv{x)]e |^ 



< -l-oo , 



^ Note that, as in a Poisson system there is no correlation of the 
particle positions, the PDF of the total field at a spatial point is 
the same whether this point is assumed to be occupied or not. 
This is a feature specific to Poisson systems, which greatly simply 
the calculation of the force PDF compared to other cases (see, 
e.g., discussion in [Ssll ). 



the above limit of Eq. (|TT|) yields 



dq 
2^ 



where we have called P{Ff^) the asymptotic shape of 
P{F^- N) for N, L,L' ^oo with N/{L + L') = no and 

g{q;no,g,fi)^~2no da; [l - cos (gge-^^)] (12) 
Jo 

is the cumulant generating function of the stochastic force 
F^. The cumulant of /'^ order of F^ is given by 



Ci{no,g,n) = i 



.1 d'^g{q;no,g,fi) 



dq' 



for ? > 1 



(13) 



g=0 



Note that Ci = (f^) and C2 = (F^) - (F^)^ and that 
for a Gaussian variable Q is & quadratic function of q. In 
our case we have that all the odd cumulants vanish by 
symmetry, while the even ones are 



C2/(?^o,g,M) 



"off 
III 



(14) 



IV. RENORMALIZATION SCHEMES 

Let us now analyze the behavior for — > 0, i.e., when 
the screening length diverges. This is equivalent to apply- 
ing the Jeans' swindle in the usual way. First of all we see 
that, as anticipated in the previous section, all the non- 
zero cumulants diverge (^ 1/m) ^-nd thus the stochastic 
force F becomes ill defined in the following sense: its 
PDF P{F) at smaller and smaller /i becomes broader and 
broader vanishing point-wise at all F. This is analogous 
to the behavior of the PDF of the sum of N identically 
distributed and independent random variables with finite 
variance and zero mean: it becomes broader and broader 
as N increases and vanishes point-wise at all finite values 
for — > 00. We underline, that this limit /i — > is also 
analogous to the one behind Chandrasekhar's calculation 
[2l| in d = 3 leading to the Holtzmark distribution; and 
likewise that considered in obtaining a well defined PDF 
for the ID case for a class of perturbed lattices in [20j . 

It is clear from Eq. (ITil) that it is possible, however, 
to converge to a well defined PDF in the limit /i — 
if one appropriately renormalizes also the characteristic 
scales of the system. There are essentially two different 
possibilities: 

• no, /i — > with uq/ fi — a > and fixed g > 0: in 
this case all the cumulants, and therefore the PDF 
P{F) are well defined with 
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21 



I 



> yi>i. 



This limit is non Gaussian as even cumulants of or- 
der larger than two are strictly positive. We call 
this case the sparse limit (SL) as the particle den- 
sity vanishes together with /i. 
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• (7,/i — with / ^JL = 6 > and fixed n-o > 0: in 
this case also P{F) is well defined. However now 
F becomes a Gaussian variable for which 

6*2 = 6y and Ca/ = for / > 2 . 

In other words in this limit we have 
1 



P{F) = 



-FV(6no) 



This limit corresponds to the limit in which a sort 
of central limit theorem applies to the force F. We 
will refer to it as the weak interaction limit (WIL). 
This is the limit which is considered in [291 . 



Thus while the "bare" stochastic force F as defined by 
the simple sum ([5]) is ill defined in the infinite volume 
infinite force range limit, the two "renormalized" forces 
(f) = Fy^ h/uq and Tp = F^/JI/g are well defined stochastic 
variables in the same limit. However only the latter is 
Gaussian. 

The first limit is manifestly quite trivial: in the system 
we consider there are just two characteristic length scales 
Ai = I/uq (typical interparticle distance) and A2 = 
(range of the interaction). Therefore a system in which 
these two lengths are substituted by two proportional 
ones X[ = kXi and A2 = fcAa is just a spatial rescaling 
of the previous system, which consequently has the same 
force F statistics P{F). As the latter is non-Gaussian for 
finite /i it is likewise so in the asymptotic limit no , /i — >■ 
with fixed ratio no//i. For the other limit, yit, (7 — > with 
fixed 5^//Lt, this is not the case: rescaling /i and g to 
fi' — kjjL and g' = \fkg does not produce a simple spatial 
rescaling, i.e., a physically equivalent system. 

To understand better the meaning of these different 
limits it proves instructive to study the spatial correlation 
properties of the force. This is the subject of the next 
section. 



V. FIELD-FIELD CORRELATIONS 

Let us consider now, rather than the PDF, the corre- 
lation properties of the force field at two distinct spatial 
points. It is instructive to do so first for the original un- 
regularized ID interaction in the finite Poisson system of 
N particles randomly distributed in the interval \—L, L'\. 
Using Eq. ([5]), we can write 

{F{x)F{y)) ^ g^ sgn{x^ - x) sgn{xj - y)j . 

Using Eq. ^ it is straightforward to show that 



{F{x)F{y)) 



N 



U + L 



g'[L' + L-2\x-y\_ 



(15) 



■ ^^if^^2 s\L' - Lf - 2{x + y){L' - L) + Axy] 



From this formula it is clear that, as for the one- 
point properties of the unregularized force, this two-point 
quantity is ill defined in the limit N,L,L' 00 with 
N/{L' + L) — uq. Note that taking the limit symmetri- 
cally (i.e. with L' — L) removes the quadratic and linear 
divergences in the second term, but leaves a linear diver- 
gence in the first term. Further, in this case, all spatially 
dependent terms are in fact finite. This is a crucial point 
in the discussion of the renormalisation below, which we 
now formulate for convenience using the smooth regular- 
ization procedure used above. 

Using instead the pair interaction defined by Eq ([7]), 
it is straightforward to obtain, using Eq. ([3]) , and taking 
the hmit N, L, L' ^ 00 with N/ {L' + L) — n^, the result 
that 



+ G^{x~y)^{F^{x)F,{y))=nog^{^i-'-\x-y\)e-^^\--y\ 

(16) 



which, as one would expect, depends only on |a; — ?/|. 
Thus the Jeans' swindle, given by the /i — >■ limit, leaves 
the expression ill-defined. Like Eq. (fT5|) with L = L' 
the term which diverges is space-independent, while the 
space-dependent part is finite in the same limit. 

It is instructive to give the results for the two-point 
properties also in reciprocal space. The Fourier trans- 
form (FT) with respect to u = {x — y) of the two-point 
correlation function G^{u) corresponds to the power 
spectrum (PS) of the total field, which is thus given by 



(fc2 + ;i2)2 



(17) 



This result can alternatively be obtained by calculating 
directly the FT of F^{x), as in Appendix El The integral 
of the PS of a stochastic field is equal, by definition, to its 
one point variance, i.e. the integral over k of Eq. (jl7p is 
equal to Gp(0) = C2(no, (7,/i). The divergence of G^(0) 
as /X — > corresponds in fc-space to the non-integrability 
of the PS (cx l/fc^) at small fc, i.e., due to the ill de- 
fined large distance correlation behaviour We note 
that the result Eq. pT]) for the PS of the screened grav- 
itational force field is in fact valid (up to a constant) 
in any spatial dimension. In three dimensions, however, 
the unscreened limit ~ l/k"^ is integrable at small k (but 
non-integrable at large k). This is an equivalent way of 
explaining why the Jeans' swindle (formulated using an 



* We note that in any dimension, for the screened interaction, 
at small k <^ /i (i.e., at scales much larger than the screening 
length the PS Sp(k; fi) ~ fc^. This means that F{x) is a 

so-called superhomogeneous |34|| . or hyperuniform 3^, stochastic 
field. The main properties of this class of stochastic fields are that 
f—^ Gfj,{x) = 0, and that the fluctuations of the field are sub- 
poissoniann, i.e., in one dimension the normalized fluctuations in 
a region of size I decrease more rapidly than For the case 

of a PS proportional to fc^ they decay as ^ which is the 

most rapid possible decay for any proper stochastic process in 
one dimension 
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exponential screening) does not work in one dimension 
while it does in three dimensions. Note, however, that 
the divergence of the variance alone does not imply in 
itself that the PDF itself is undefined, which we have 
shown to be the case in one dimension. Indeed the total 
force PDF in three dimensions (i.e. the Holtzmark distri- 
bution) has infinite variance, due however to the singular 
behaviour of the 3D pair interaction at vanishing separa- 
tion and not to the large scale contributions. 

Let us now consider again the renormalization schemes 
introduced in Sect. IIVI It is clear that in both cases we 
obtain 



(18) 



or, equivalently, 



SF{k;fi) 2TTC2S{k), 

where stands for the (finite) field variance after 

the renormalization procedures. These procedures thus 
keep finite the dominating and diverging contributions to 
Eq. ([15]), but at the same time send to zero all other sub- 
dominant contributions. In particular this means that 
they send to zero all terms depending on the spatial ar- 
gument of the correlation function. In other words, while 
the regularization-|-renormalization procedure, in both 
schemes, makes F well defined as a one-point stochastic 
quantity, the by-product is to eliminate any space varia- 
tion of this field. The field which remains is a finite ac- 
celeration off-set of the whole system. From a dynamical 
point of view this corresponds simply to a translation of 
the whole system, and these renormalisation procedures 
thus erase all information about the relative motion of 
particles. Indeed, in the spirit of Mach's principle, all 
relevant physical information about the force field is lost 
in this limit: "A particle's inertia is due to some (unfor- 
tunately unspecified) interaction of that particle with all 
the other masses in the universe; the local standards of 
nonacceleration are determined by some average of the 
motions of all the masses in the universe, [and] all that 
matters in mechanics is the relative motion of all the 
masses" [36j . 

These observations naturally lead us to consider the 
statistics of the differences in forces between spatially 
separated points in the infinite system limit. 



VI. RELATIVE FIELD ANALYSIS 

Let us begin again by considering the homogeneous 
Poisson distribution defined by Eqs. ([T]) and ([3]) of 
particles in the interval [—L, L'] and interacting by the 
pair force field (U). Let us fix an arbitrary point, say the 
origin x — 0, and consider the difference of the total field 
at a point x and the origin: 

N 

SF{x) = F{x) - F{0) = 2gY^[0{x^ ~ x) - e{x,)] . (19) 



As the sum is simply, up to a sign depending on x, equal 
to the number of particles in the interval [0,a;), 6F{x) 
is manifestly independent of the extremes of the interval 
[—L, L'] and therefore remains the same and well defined 
in the limit N, L,L' ^ 00 with fixed N/{L+L') = no > 
taken in any arbitrary way. It is thus simple to calculate 
the one-point PDF of 6F{x), using the properties of the 
Poisson distribution: 



P{SF;x) 



-nokl KN)' s[SF + 2;g(sgn(x))] (20) 
1=0 



In other words SF{x) can take only (positive or negative 
depending on the sign of x) integer multiple values of 
2g with a Poisson probability distribution of mean no 1 2;]. 
Proceeding as for the analysis of P{F) in Sec. IIIIl using 
the identity Eq. dTU]) . we can derive the characteristic 
function P{q;x) — FTsfP{SF;x) (where FTsf indicates 
the FT with respect to SF) as 

P{q;x) = exp [Qiq;x)] 

where the cumulant generating function Q{q] x) of SF{x) 
is given by 



-no 



\x\ (1 



J2gq sgn( 



(21) 



Using the definition Eq. (jl3p we can obtain the cumulants 
Xj(x) of SF{x) for all j > 1 as 



^ji^) = [-2gsgn(a;)]^no|a;| 



(22) 



To calculate the same quantity, but now using the 
"smooth" Jeans' swindle formulated as the /i — > limit 
of the screened gravitational interaction Eq. ([7]), we fol- 
low, as in Sec, mil the Chandrasekhar derivation starting 
directly from the definition of the PDF: 



P,.{5F;x) 



N 



-L' 



dxi 
L + L' 



SiSF 



N 
i=l 



n 

[sgn{xi - a:;)e"^l^"~^l - sgn(a;i)e~'^l^' 



(23) 



Following the analogous manipulations and taking the 
limit N,L,L' 00 with fixed N/{L + i') = "0 > 
as in Sec. IIIIl we obtain that the cumulant generating 
functional can be written as 



-no 



dy X 



sgn(a)e 



and thus the cumulants are given as 
Ki^) = i-gyno / dy X 



(24) 

'0 

(25) 



sgn(?/)e — sgn(?/ — x)e^ 
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First of all we notice that 

A5' (x) =0 Va* > 

It is, furthermore, straightforward to verify that for inte- 
ger I 

while 

Finally one can show that for j > 2 

limX^{x)^X,{x) (26) 

i.e., the cumulants of order j > 2 converge, for /i — > 0+, 
to those derived above in Eq. ((22|) for the case without 
screening, while the average value (j — 1) instead van- 
ishes: the Jeans' swindle simply removes the average den- 
sity, thus making the average force zero everywhere. The 
associated generating functional is thus 

Qiq; x) = -no\x\ (l -I- 2igqsgn{x) - e'^aqs^nix)^ (27) 

Note that this result may be obtained directly from 
Eq. (PO)) by simply replacing £ inside the delta function 
by ^ — no I x I, i.e., by simply subtracting by hand the con- 
tribution of the mean density no to the difference in the 
force SF{x). Note that, conversely, we can also obtain 
the initial result (without the Jeans' swindle) using the 
second derivation, but putting ^ = in Eq. (|25|) before 
doing the integral. In other words, the two limits, (i) the 
extremes in the integral in Eq. (|25|) and (ii) — >■ 0, can- 
not be exchanged. The price to pay for this exchange is 
the uniform contribution coming from the mean density. 

It is straightforward also to calculate the two-point cor- 
relation functions of SF{x). By using again the joint par- 
ticle positions PDF of Eq. we can simply evaluate the 
averages and then take the thermodynamic limit which 
is now well defined, finding 

Gdiff(a;,y) = (<5F(a;)5F(y)) = (28) 
2g^no[\x\ + |y| - \x - y\ + 2no\x\\y\] 

for the unregularized case (i.e. without Jeans' swindle), 
and 

Gdiff(x,y) = 2g^no[\x\ + \y\ - \x - y\] (29) 

when the Jeans' swindle is used. The latter result is most 
easily recovered by calculating the correlation function at 
finite fi using the result in Eq. (|16|) . and then taking the 
limit /i — ?> 0. The additional quadratic term in Eq. ([25)1 
is simply the contribution from the non-zero mean den- 
sity. The interpretation of the other terms, common to 
both expressions, is very simple: (i) when x and y have 
different signs, the two intervals [0, a;) and [y, 0) have 



empty intersection, and, as there is no correlation be- 
tween the position of particles in a Poisson distribution, 
the fluctuations in the variables SF{x) and SF{y) are sta- 
tistically independent, and the correlation function [the 
"connected" part for Eq. (|28p ] is therefore zero; (ii) if in- 
stead X and y have the same sign, the segments [0, x) and 
[0, y) overlap with intersection equal to the shorter of the 
two segments, and thus there is an non-zero correlation, 
proportional to the length of this interval. 



VII. CONCLUSIONS 

In conclusion let us consider the implications of our 
results for the question of whether an infinite system limit 
may be defined for the dynamics of a ID system of points 
interacting by the ID version of Newtonian gravity, when 
the initial distribution of these points is Poissonian. More 
specifically we wish to consider the dynamics of particles 
described by equations of motion in one dimension given 

X, + 2Hxi = -4 1™ H ^Sn(a;z - a;j)e~^l^'-^^l, (30) 

where a{t) is a function describing the expansion of a 3D 
universe, and H = d{\n a) / dt is the corresponding Hubble 
expansion rate (and the case a{t) = 1 describes the static 
universe limit). While in the limit ^ ^ {) these equations 
are explicitly well defined for a finite number of particles, 
the question is whether they remain well defined when 
we consider the usual thermodynamic limit (as defined 
above: L, L' — >■ cx) at constant no = N/L -\- L'). The 
importance of this limit is that it models the case of an 
infinite universe, which is the application of relevance for 
these toy models. 

We have shown in [l^] that this limit may indeed be 
defined in the case of an infinite array of particles ini- 
tially displaced off a perfect lattice, for a broad class of 
such displacements. More specifically we did so by calcu- 
lating the PDF of the force as defined on the right-hand 
side of Eq. (|30l) . and showing it to be well defined for 
this class of distributions. In this paper we have instead 
shown that the analogous PDF is not well defined for the 
case of a homogeneous Poisson particle distribution in 
the same limit, and thus that the infinite system limit is 
not defined in the same sense. 

In our discussion of other possible regularisations of 
this limit, we have shown that the divergence in the total 
force arises from the dominant contribution of particles 
infinitely far away. Because the ID gravitational force 
is independent of separation, however, this component 
does not contribute to the relative force on any two par- 
ticles at a finite distance. As a result, while the force at 
any point itself becomes completely undetermined, the 
difference between forces at two spatially separate points 
does not. This means that while Eq. ([50)1 is a badly 
defined equation of motion for each particle i, one can 
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nevertheless write a well defined equation for the relative 
displacements Sxij = Xi — Xj of two particles i and j: 

6x,j + 2HSx,j = lim [F^{x,) - (31) 

where 

F^{x^) = -4 E «g^(^^ - 2:fc)e-''l^'-^''l . (32) 

Thus, if we consider the evolution from homogeneous 
Poissonian initial conditions, the position of a particle 
after any finite time will always depend on /i, and di- 
verge with probability 1 as /i — > 0. On the other hand, 
the relative position of any two particles initially at a 
finite distance will extrapolate to a finite //-independent 
value in the same limit. In other words, the clustering 
dynamics of the system — entirely characterised by the 
relative positions of particles (e.g. by two or higher point 
connected correlation properties of the density field) — is 
well defined. In the spirit of Mach's principle, the diverg- 
ing absolute displacement, in an infinite system which has 
intrinsically no centre or preferred point, is not, in any 
case, of physical relevance. 

We have shown the above statements to be true strictly 
only at the initial time, i.e., the clustering dynamics is 
well defined in an infinite Poisson particle distribution as 
we have described. That they remain true as the system 
evolves away from the initial Poisson distribution can be 
most easily verified by considering the evolution of the 
density perturbations in fc-space. It is well known (and 
straightforward to show — see e.g. [13]) that, in any 
spatial dimension, the evolution of a self-gravitating sys- 
tem leads to a /c-independent amplification of the PS of 
density fluctuations at small k, provided the PS does not 
vanish faster than fc'' at fc = 0. The small k (i.e. large 
scales) scaling behaviour of the PS of a Poisson parti- 
cle distribution is thus unchanged by evolution. As seen 
above [cf. discussion around Eq. (IT7)) ]. it is this behaviour 
which determines the convergence properties of the force. 
These properties therefore remain invariant under evolu- 
tion, guaranteeing that the clustering dynamics remains 
well defined. 

The meaning of this limit is that it corresponds to the 
clustering dynamics of scales much smaller that 1/fi for 
the screened pair interaction, and the limit /i — > sends 
this upper cut-off scale to infinity. We remark that this 
limit exists because the time scale for evolution of clus- 
tering at a given spatial scale a; > increases with x, 
i.e. the clustering is what is known as "hierarchical", 
proceeding from the smallest scale upwards. It is easy 
to show heuristically, as follows, that this corresponds, 
in general, to a condition on the scale dependence of the 
relative forces which is indeed satisfied here. The charac- 
teristic time scale t^ for evolution of a system on a scale 
X can be estimated as 

^ \IW{xj\ ^^^^ 



where |AF(2;)| is the typical relative force on points at 
separation x. For the case of ID gravity in an infinite 
Poisson distribution, employing the Jeans' swindle as 
above, we have {[F{x) — F{0)]) — and therefore we 
take |AF(x)|2 = {[F{x) - F{0)]^) = Gdiff (2;, a;). From 
Eq. (|29p it follows then that t^ ~ x^^^. Note, however, 
that if we do not employ the Jeans swindle, we have 
instead that the typical force difference on a scale x is 
I {[F{x) - F(0)]) I = 25^00; > 0. Therefore now t^ is in- 
dependent of scale, i.e., all scales evolve on the same time- 
scale, which means that the infinite system limit cannot 
be defined. This is, indeed, the fundamental physical 
reason for the introduction of the Jeans' swindle: it re- 
moves the centre of the system towards which all points 
otherwise collapse in a time independent of the system 
size. 



Finally let us comment on related results given recently 
in j29|l. This paper derives, using the Jeans' swindle in 
its usual formulation as a symmetric "top-hat" sum, the 
PDF of the gravitational force in one dimension, defining 
it taking the equivalent of the rcnormalised weak coupling 
limit we have discussed. A central point in the paper is 
the observation that in passing from d = 3 to c? = 1 the 
statistics of the gravitational force as characterised by 
the PDF changes from the power-law tailed Holtzmark 
distribution to the Gaussian form obtained in one dimen- 
sion Given our results and discussion here, we consider 
that there is no basis for giving any significance to this 
fact: the Gaussian PDF in one dimension is not obtained 
in the analagous limit to that used in three dimensions, 
and the modified renormalized limit which gives it has 
only a trivial physical significance as it leads to a spa- 
tially trivial force field. Further we note [2l|, [2^ that the 
"fat" (i.e. non-integrable) power-law tail of the Holtz- 
mark distribution in fact arises from the divergence of 
the pair interaction at small separations, and has nothing 
to do with its long range nature. Indeed, even without 
regulation of the singularity at small separations, other 
distributions in three dimensions (e.g. "shuffled" lattices 
with exclusion regions around each particle (2 61] ) have a 
Gaussian tail in the gravitational force PDF. More gen- 
erally, in fact, as we will discuss in forthcoming work [38{. 
the long-range nature of a pair interaction does not lead 
to divergences of all moments of the force PDF of order 
larger than a typical value, leaving the PDF itself defined 
even though power law tailed. This is exemplified in the 
case we have analysed: from Eq. (ITU) we see that the 
cumulants of the force field diverge at any order at the 
same rate as /i 0. 



We acknowledge Frangois Sicard for many fruitful dis- 
cussions in the context of our collaboration on related 
projects. We also thank Bruno Marcos and Bruce Miller 
for useful conversations. 
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Appendix A: Derivation of the power spectrum of 
the field F{x) 

The power spectrum Spik; /i) given in Eq. (|17|) can be 
also derived directly from its definition: 



lim 

L,L'— >oo 



\F^ik;L,L')\' 
L + U 



(Al) 



where 



L' 



with F^{x) given by Eqs. (I7][8]) and k an integer multiple 
of 27r/(L + L'). In order to evaluate F^{k\L,L') in a 
simple way it is useful to notice that the pair force ([7]) 



can be derived by a simple derivative (and a change of 
sign) from the pair potential 
is a solution of 



ge '^I'^l / ^ which 



dx^ 



(A2) 



The FT (in [~L,L']) ^{k;L,L') of the total poten- 
tial associated with a density field n[x) thus satisfies 
{k'^ + fj.^)^{k;L,L')) = -2gn{k] L, L'). Now substituting 
\F{k;L,L')\^ = k'^\^{k;L,L')\'^ in Eq. and using 

the definition of the PS of a Poisson point process 



|n(fc;L,i')r 
lim — = no , 

L,L'^oo L + L' 



we obtain Eq. p7)) . 
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